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Abstract — Charge migration through DNA is a problem 
to be solved. A numerical wavefront solution for quantum 
transmission lines with charge discreteness is obtained 
as model for the charge migration of a chain DNA 
molecules induced by electromagnetic radiation. The 
nonlinearity of the system becomes deeply related to 
charge discreteness. The wavefront velocity depends on 
the normalized (pseudo) flux variable. Finding the 

dispersion relation for the normalized flux (|) / (|) 0 we show 

that the condition V I. 2 > 0 on the wave -front velocity 
gives the band- gap conditions for the charge propagation 
on the system.. 
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I. INTRODUCTION 

Currently nanostructures are embedded in countless 
devices and systems [1-5] . Naturally, at this scale and for 
low temperature, quantum mechanics plays a key role. 
Lately, much effort has been dedicated to study 
nanostructures, using as model of quantum circuits with 
charge discreteness [6-12]. On the other hand, the 
discovery of charge migration in deoxyribonucleic acid 
(DNA) stimulated intensive investigations of the 
electronic properties of DNA due to their significance in 
biosynthesis and radiation-induced damage and repair 
processes [1-3]. Furthermore, considerable interest in 
nanodimensional structures of DNA possessing unique 
self-assembling and self-recognition properties has 
increased the last decade in connection with the 
possibility of the development of molecular 
nanoelectronic devices which are expected to provide 
high storage of information and high-speed signal 
processing within a wide temperature range [4-6]. In fact, 
DNA molecules can be well combined with silicon 
technology transcending the potential of the present 
quantum wires and are supposed to be used in modern 
computer technology as a binary data structure by 
applying a programmable linear self assembly of the 
sequence of complementary nucleic base pairs of DNA. It 


is in this context that we are interested in spatially 
transmission lines with charge discreteness [13]. In this 
work we will consider a wavefront solution for a quantum 
circuit (transmission line) simulating the DNA molecule 
in a cell sample. A numerical solution is founded and 
characterized for this specific system with charge 
discreteness, which can extended for the description of 
more complex extended systems a. In section 2 we will 
introduce a generalization of classical transmission line 
with charge discreteness and we present the quantum 
transmission lines with charge discreteness, the 
Hamiltonian for coupled circuits, and the equations of 
motion for the spatially continuous system. In Sec. 3, the 
wavefront solution is considered. Finally, we give our 
conclusions. 

II. GENERALIZATION OF TRANSMISSION 
LINE AND QUANTUM DISCRETENESS 

For a chain of molecules of DNA, we consider a 
homogeneous classical transmission line, assumed 
infinite, where the every cell is constituted of a LC circuit 
with inductance L and capacitance C per unit length. 
Assume that the interaction between neighbor cells is 
trough the capacitors (direct line), the classical evolution 
equations for the electrical current and charge, become in 
this case 

L ^ im= d (2qm_qm+1_qm_l) (1) 

Where the integer m designates the cell at position in the 

d 

chain and, as usual, 1 = — q m is the electrical current. 

’ ’ m t , T.m 

dt 

The above linear equation of evolution becomes directly 
from the classical Lagrangian L ag given by 

Using the Lagrange equations one obtains (1). 

Supporting solutions like a plane wave. Explicitly we 
have, 
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q in =q 0 exp(icot-ikm) 


(3) 


Where, as usual in this case, the frequency 00 and the 
wavenumber k become related through the dispersion 
relation. The variable q 0 is a constant parameter 

describing the amplitude of the wave. Note that we have 
considered by sake of simplicity the lattice-size a as one, 
(a =1). That is the adimensional number kin the 
equation (6) is really ka . 

After some algebra, one obtains from (3) the dispersion 
relation. The infrared limit ( ka — > 0 ) related to the 
continous -space structure becomes direct and 


corresponds to the dispersion relation 
namely a non dispersive medium. 


© = k/ VLC , 


From CO = 
velocity 


k/ VLC , we get the phase velocity and group 


v =“=i/7Ec,v =^=i/Vlc (4) 

p k s dk 

From a general point of view, for arbitrary composition of 
a cell in the line, the Hamiltonian of this generalized 
electrical transmission line becomes quadratic, namely, 

H =z4 < f ) "-- < I’"A (5) 

m s Z Z U 


As illustration and following [6-12, 14], and from the 
Hamiltonian (5), the usual quantization procedure for flux 
and charge, and the prescription (2) for charge 
discreteness, we could construct the quantum Hamiltonian 
for the direct transmission line with charge discreteness ( 

q e ), which may be written as: 




2r . 2 q 1 2 

u sm n 


where the index m describes the cell (circuit) at position 
m, containing an inductance L and capacitance C. The 


conjugate operators, charge q and pseudoflux (]) , satisfy 
the usual commutation rule 


Qm’ 


= i ^ 5 m,m and 


q m >q s 


•LA 


= o. 


(7) 


A spatially extended solution of Eq. (6) corresponds to 
the quantization of the classical electric transmission line 

with discrete charge (i.e. elementary charge q e ). Note that 


in the formal limit q e — » 0 the above Hamiltonian gives 

the well-known dynamics related to the one -band 
quantum transmission line, similar to the phonon case. 
The system described by Eq. (10) is very cumbersome 


since the equations of motion for the operators are highly 
nonlinear due to charge discreteness. However, this 

system is invariant under transformation q k — » q k + q , 

that is, the total pseudo flux operator (j) = E*. 

commutes with the Hamiltonian; simplifying the study of 
this system. 

H m = sm ^<|> m - — (q m+1 -<|> m ) (8) 

Lq e zn 2C 


where H m represents the Hamiltonian density operator 

for the fields. From the above Hamiltonian we find the 
equations of motion (Heisenberg equations) for the field 
operators: 

0 1 

^m=-(qm + l+qm-l- 2 q m ) (9) 

■— q = — — — sin(^- 4>) (10) 

ot Lq e h 

and from eq. (6), the dispersion relation is given by 
|sink/2| 

co = 2 J — (11) 
VLC 

If kD 1, © = k/ VLC , from here we can obtain 


v =^ = l/VLC,v =^ = 1/^C(12) 

p k g ak 


III. WAVEFRONT SOLUTION 

In real seismic applications there is always the presence 
of damping. We shall consider the effect of its simplest 
form, small viscous damping. Eq. (3) is extended by 
adding a linear damping term X : 

Now we consider a homogeneous classical transmission 
line, assumed infinite, where the every cell is constituted 
of a LC circuit with series inductance L and shunt 
capacitance C per unit length . The approach is based on 
the mapping of field components (i.e., E and H) in the 
medium to the voltages and currents of the equivalent 
distributed L-C network [15, 17]. It is well known that 
dielectric properties like permittivity and permeability can 
be modeled using distributed L-C networks. To illustrate 
how these material parameters relate the distributed series 
impedances and shunt admittances of the network, 1-D 
distributed L-C network is depicted in Fig. 1. The 
network consists of a series per-unit-length impedance Z ' 
in z direction and a shunt per-unit-length admittance Y' in 
y direction. 
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Fig.l: Transmission line model with L , C per unit length. 


For the conventional homogeneous isotropic RHM with 
positive permittivity 8 and permeability p, Eqs.(l, 2) 
implies a network of a low-pass topology with series 
inductor L ' = p (H/m) and shunt capacitor C ' = s (F/m), 
or Z=jodL ' and Y = /coC ' , both of which are positive 
quantities. The negative permittivity and permeability 
[16, 17] in LHM leads to the question whether the 
parameters Z ' and Y ' in the network representation can 
also be made negative. From the impedance perspective, a 
negative L ' and C ' can be realized by exchanging their 
inductive and capacitive roles, which means the series 
inductor becomes a series capacitor, and the shunt 
capacitor becomes a shunt inductor. In this paper we are 
interested in the normal RHM transmission line which in 


the limit q e — >0 becomes the usual one 
transmission wave equation. 

V 2 = CO 2 / k 2 = 1 / LC . In general the dispersion 
relation from Eqs. (9) and (10) is 


cd 2 LC 


= cos(q e (| )/h)( 


sin(k IT) 
k/2 


) 2 


where V is the phase velocity given by 

I cd 2 LC 


V — < 


=V cos (q e (i )/ft)( sin(1 ^ 2) 

k/2 


band 


(13) 


(14) 




Fig.2:. The normalized expression (17) 


is shown as a 


function of (|) / (|)q and k 


i—i 



Fig. 3: Profile of \ 2 as a function of c|) / c()q with k = cte 


Fig. 4: Profile of\ 2 as a function of k with ([)/ c|)q = cte 

2 

Figures 2-4 show the normalized V of the charge 
migration in a transmission line where we can see that the 
velocity of wavefront which has bands and gaps. To see 

2 

this in a clearer way we know that always v must be 
positive (V 2 LC>0), so the two equations (13), (14), 
are satisfied simultaneously when 4> / 4 >q = 2n7l , 

n = ±1, ±2... that means that (|) = 2nn§ 0 = nh / q e . 
Let 

. 2t ^ ® 2 LC . , . .sin(k/2). 2 
f (k) = v 2 LC = — — = cos(pk)( v 0 2 (15) 

R R / Z 
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Where we define pk = q e (|) / tl — (|) / (|) 0 , so we have 
the factor p as a free parameter to analyze the quantum 
circuit. 

In figure 5 we present the case with (J)/ (J) () =0 where 
f (k) > 0 and there is not a structure of bands and gaps 
because q e — >0. In figure 6 we have a structure of 


bands and gaps because pk = q e (|)//2 = (|)/(|) 0 ^ 0 



Fig. 5: Plot of the velocity of the wavefront , a function 

of k with the magnetic flux as the parameter p. 


As specified by Eq. (18), if p=0, there is only a structure 
of bands because pk = q e (|) / tl — (|) / (|) 0 = 0 , that is 

q e ->o. 



Fig. 6. Plot of the velocity of the wavefront , as a function 
of k with the magnetic flux as the parameter p. 

As specified by Eq. (18) there is a structure of bands and 

gaps because pk = q e (J> / h — (|) / (|) 0 ^ 0 with p=3 , so 

that have five bands where there is wave propagation 
and charge migration as are shown in Figure 7. 



-r 


Fig. 7. V 2 LC > 0 where there are wave propagation. 
As we have seen from these graphs, the condition 

V 2 > 0 on the wave-front velocity gives the band-gap 
conditions on the system. In fact, from (16) the restriction 


LC 


cos(q e <|)/ K)( 


sin(k / 2) 
k/2 


) 2 >0 


means that there exists a sequence of bands and gaps. 

The nonlinearity of the system becomes be deeply related 
to charge discreteness or in terms Of discrete values of 

quantum flux (|) 0 = ti / q e . The wavefront velocity is 


found to depend on a step discontinuity on the (pseudo) 
flux variable, p, displaying allowed and forbidden regions 
(gaps), as a function of p. 

Defining, (|) / (|) 0 = pk , any complex function ¥ can be 


written as \\f = \|/ Q exp(ipk) where \|/ Q is the amplitude 

and 0 is the phase. It is obvious that changing the 
phase pkby 2n7Twill not change \| J and, 

correspondingly, will not change any physical properties 
in this quantum transmission line. However, in the 
superconductor of non-trivial topology, e.g. 
superconductor with the hole or superconducting 
loop/cylinder, the phase pkmay continuously change 

from some value pk Q to the value pk 0 + 2TCI1 as one 

goes around the hole/loop and comes to the same starting 
point. If this is so, then one has n magnetic flux quanta 
trapped in the hole/loop, [18, 19], In the case of charge 
migration in DNA, charge discreteness is the statement 
that charge comes in packets which are of size 1 electron 
charge. These results enable the investigation of charge 
migration as quantum charge in both longitudinal and 
transverse configurations and stimulate theoretical 
interpretations [20-24]. Some of these studies are 
concerned with the issues of the charge migration induced 
by environmental factors, among which ionizing radiation 
or RF fields are of great interest. 

Charge migration through DNA has been the focus of 
considerable interest in recent years because DNA is the 
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molecule that contains all of the information required to 
build and maintain the cell . A deeper understanding of 
the nature of charge transfer and transport along the 
double helix is important in fields as diverse as physics, 
chemistry and nanotechnology. It has also important 
implications in biology, in particular in DNA damage and 
repair. In this context this paper can be a contribution in 
the study of this topic. 

IV. CONCLUSIONS 

A charge migration in a chain of DNA molecules is 
simulated with the quantum electric transmission line 
with charge discreteness described by the Hamiltonian 

(12), and equations of motion (13-14), Wavefront 
solution was found.. One condition on the velocity 
generates a band-gap structure dependent on the pseudo 
flux parameter, namely, there exist regions for which a 
solitary wavefront propagates with constant speed 

according the value of (j) / (j)^ . The main results of this 
work are the existence of the band-gap structures for 
0 < (j) / < 1 . Charge migration under radiation is 

possible when v 2 > 0 . 

REFERENCES 

[1] S. Datta, Lessons from Nanoelectronics, a Lecture 
Notes Series (Volume 1, Publication, Singapore 
: World Scientific Press, 2012). 

[2] S. Lyshevski, Nano and Molecular Electronics 
Handbook, CRC Press. Taylor & Francis Group, 
2016. 

[3] H. J. de Los Santos, Principles and Applications of 
Nano MEMS Physics. Springer, Berlin, 2005. 

[4] David K. Ferry and Stephen M. Goodnick, Transport 
in Nanostructures. (Cambridge University Press, 
2009). 

[5] W. D. Heiss (Ed.), Quantum Dots: Doorway to 
Nanoscales Physics. (Springer, Berlin 2005). 

[6] Y. Q. Li and B. Chen, Phys. Rev. B 53, 4027 (1996). 

[7] B. Chen et al, Phys. Lett. A 335, 103 (2005), and 
references therein. 

[8] J. C. Flores, Phys. Rev. B 6 , 235309 (2001). 

[9] J. C. Flores, Europhys. Lett. 69, 116 (2005). 

[10] T. Lu and Y. Q. Li, Mod. Phys. Lett. B, 16, 975 

( 2002 ). 

[11] S. A. Gurvitz, Phys.Rev.Lett.85, 812 (2000). 

[12] C. A. Contreras and D. Laroze, Physica B, 476, 77- 
81, (2015). 

[13] N. V. Grib et al, Distance-dependent coherent 
charge transport in DNA: crossover from tunneling 
to free propagation, Journal of Biophysical 
Chemistry, Vol.l, No.2, 77-85 (2010). 


ISSN: 2349-6495(9) / 2456-1908(0) 

[14] J. C. Flores et al, Phys. Rev. B 74, 193319, (2006), 
and references therein. 

[15] H. Torres-Silva, "A Transmission Line Model with 
Metamaterial Effects in Gamma Ray Bursts", 
Journal of Electronics and Communication 
Engineering Research, Volume 3, Issue 5, P 12-17, 
(2016). 

[16] H. Torres-Silva, D. Torres-Cabezas and J. Lopez- 
Bonilla, Charge discreteness in LC quantum circuit, 
Open Journal of Technology & Engineering 
Disciplines (OJTED), Vol. 2, No. 3, September 

2016, pp. 01-06. 

[17] H. Torres-Silva, Metamateriales, Actualidad y 
Desarrollo, Dyna, vol 92, n° 3, Cod 8244, 2017, and 
references therein. 

[18] S. Frankel, National Mathematics Magazine, Vol. 
11, No. 4, pp. 177-182(1937) 

[19] D. Bascom; F. William , Review Letters. 7 (2): 43- 
46, (1961) 

[20] R. Doll, N. Nabaue, (July 1961).. Physical Review 
Letters. 7 (2): 51-52, (1961). 

[21] B. Giese B, S. Wessely, M. Spormann, U. 
Lindemann, E. Meggers and M.E. Michel-Beyerle, 
On the mechanism of long-range electron transfer 
through DNA, Angewandte Chemie International 
Edition 38(1999), 996-998. 

[22] V.D. Lakhno, Soliton-like solutions and electron 
transfer in DNA, Journal of Biological Physics 26 
(2000), 133-147. 

[23] S. Eh. Shirmovsky and D.L. Boyda, Study of DNA 
conducting properties: Reversible and irreversible 
evolution, Biophysical Chemistry 180-181 (2013), 
95-101. 

[24] C.T. Shih, Sequence and energy dependence of 
electric transport properties of DNA -A tight- 
binding model study, Chinese Journal of Physics 45 
(2007), 703-707. 


www.ijaers.com 


Page | 211 


